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ABSTRACT 

The discovery of very high energy (VHE) gamma-ray emitting X-ray binaries has triggered an intense effort 
to better understand the particle acceleration, absorption, and emission mechanisms in compact binary systems, 
which provide variable conditions along eccentric orbits. Despite this, the nature of some of these systems, and 
of the accelerated particles producing the VHE emission, is unclear To answer some of these open questions, 
we conducted a multiwavelength campaign of the VHE gamma-ray emitting X-ray binary LS I -1-61 303 in- 
cluding the MAGIC telescope, XMM-Newton, and Swift during 60% of an orbit in 2007 September. We detect 
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a simultaneous outburst at X-ray and VHE bands, with the peak at phase 0.62 and a similar shape at both 
wavelengths. A linear fit to the simultaneous X-ray/VHE pairs obtained during the outburst yields a correlation 
coefficient of r = 0.97, while a linear fit to all simultaneous pairs provides r = 0.81. Since a variable absorption 
of the VHE emission towards the observer is not expected for the data reported here, the correlation found 
indicates a simultaneity in the emission processes. Assuming that they are dominated by a single particle pop- 
ulation, either hadronic or leptonic, the X-ray/VHE flux ratio favors leptonic models. This fact, together with 
the detected photon indices, suggests that in LS I H-61 303 the X-rays are the result of synchrotron radiation of 
the same electrons that produce VHE emission as a result of inverse Compton scattering of stellar photons. 
Subject headings: binaries: general — gamma rays: observations — stars: emission-line. Be — stars: individ- 
ual (LS I H-61 303) — X-rays: binaries — X-rays: individual (LS I H-61 303) 



1. INTRODUCTION 

LS I +61 303 is one of the few X-ray binaries that have 
been detected in ve ry high energy (VHE) gamma rays (see, 
e.g., iParedesI 120081 for a recent review). It is a high-mass 
X-ray binary system loc ated at a distance of 2.0±0.2 kpc 
(iFrail & Hjellmin jll991l) . The system contains a rapidly ro- 
tating early type BO Ve star with a persistent equatorial decre- 
tion disk and mass loss, and a compact object with a mass 
between 1 and 4 M(7) orbiting it every ^26.5 d in an eccen- 
tric orbit (see ICasares et all 120051: iGrundstrom et all 120071: 
lAragona et al.ll20'09[ and references therein). LS I +61 303 
was classified as a microquasar based on structures detected 
from five to several tens of milliarcseconds with the Euro- 
pean VLBI Network (EVN) and MERLIN (Massi et al. 2004, 
and references therein), although analysis of later EVN and 
MER LIN data sets does no t reveal the presence of such struc- 
tures dPhawan et all 120061: lAlbert et al.ll2008al) . Very Long 
Baseline Array (VLBA) images with a resolution of 1 mas 
(2 AU at the source distance) obtained during a full orbital 
cycle show an elongated structure that rotates as a function of 
the orbital phase (Dhawan et al. 2006). Later VLBA images 
show repeating structures at the same orbital phases, reinforc- 
ing the idea th at the milliarcsecon d structure depends on the 
orbital phase (lAlbert et al.l l2008al) . This may be consistent 
with a model based on the interaction between the relativistic 
wind of a young non-accreting pulsar and the wind/decretion 
disk of the stellar companion (Dubus 2006a), as occurs in 
PSR B 1259 -63 (but confining the pulsar wind may be prob- 
lematic; see iBogovalov et"ani2008l) . 

LS I +61 303 shows periodic non-thermal radio outbursts 
on average every Porb=26.4960±0.0028 d, with the peak of 
the radio emission shifting progressively from phase 0.45 to 
0.95, using ro=JD 2,443,366.775, with a modulation period of 
1667±8 d (Gregorv 2002). According to the most precise or- 
bital parameters, the periastron takes place at p hase 0.275 and 
the ec centricity of the orbit is 0.537 ± 0.034 (lAragona et al.l 

llool . 

LS I + 61 303 has been observed several times in X- 
rays (see ISmith et aTl 120091 and references therein). It 
generally displays quasi-periodic X-ray outbursts, with 
the m aximum occurring be t ween o rb ital phases 0. 4 and 
0.8 (Goldoni & Mereshettil Il995[ iTavlor et al." 119961: 
l^redes et al. 1997; Harriso n et al.l [2000; Esposit o et al.l 
12007), although the lack of a sensitive long-term monitoring 
has prevented to search for a super-orbital modulation. The 
source also shows shor t-term variability on timescales of 
hours (ISidoU et al.ll2006l) . 

At VHE gamma rays, LS I +61 303 has been clearly 
detected several times both by MAGIC ("Albert et al.l 120061 
l2008a, 2009) and VERITAS ( Acciari et al. 2008, 2009!). The 
source also displays VHE gamma-ray periodicity, with min- 



ima taking place near periastron, where only upper limits were 
found in MAGIC observations, and maxima occurring on av- 
erage at phase 0.6-0.7, although the source has also shown a 
secon d peak at phase 0.84 in one of the cycles (lAlbert et al.l 
l2009l) . There are indications of correlated X-rayA'^HE emis- 
sion, based on non-simu ltaneous data t aken more t han 6 hr 
apart (lAlbert et al.ll2009l) . and 1 d apart (lAlbert et al.l l2008a). 

The lack of a systematic behavior from cycle to cycle at 
X-ray and VHE bands, and the occurrence of short-term vari- 
ability in the X-ray flux, has hampered the definitive detec- 
tion of an X-ray/VHE correlation from the comparison of 
non-simultaneous data. We therefore conducted a multiwave- 
length campaign in 2007 September, covering the epoch of 
maximum VHE emission. Here, we report the first simulta- 
neous VHE and X-ray observations of LS 1+61 303 obtained 
with the MAGIC Cherenkov telescope and the XMM-Newton 
and Swift X-ray satellites that show correlated emission in 
both energy bands. We also briefly comment on radio and 
optical spectroscopic observations obtained during the cam- 
paign. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. VHE Gamma Rays 

The VHE gamma-ray observations were performed from 
2007 September 4 to 21 using the MAGIC telescope on the 
Canary Island of La Palma (28.75°N, 17.86°W, 2225 m a.s.l.), 
from where LS I +61 303 is observable at zenith angles above 
32°. The essential parameters of MAGIC are a 17 m diame- 
ter segmented mirror of parabolic shape (currently the largest 
dish for an air Cherenkov telescope), an f/D of 1.05, and an 
hexagonally shaped camera of 576 hemispherical photo mul- 
tiplier tubes with a field of view of 3.5° diameter MAGIC 
can detect gamma rays from 60 GeV to several TeV, and with 
a special setup the trigger threshold can be reduced to 25 GeV 
(Aliu et al. 2008). Its energy resolution is AE /E = 20% above 
energies of 200 GeV. The current sensitivity is 1.6% of the 
Crab Nebula flux for a Scr detection in 50 hr of on-source 
time. The improvement compared to the previo us sensitivity 
was achieved by installing new 2 GHz FADCs (lAlbert et al.l 
2008b). 

The total observation time was 58.8 hr, including 39.6 hr 
in dark conditions a nd 19.2 hr under moonlight or twilight 
dBritzger et al.l2009l) . The range of zenith angles for these ob- 
servations was [32°, 50°], with 96% of the data having zenith 
angles below 43°. The final effective total observation time 
is 54.2 hr. Table [T] gives the MJD, the effective observation 
time, and the orbital phase for each of the MAGIC obser- 
yations. The obs ervations were carried out in wobble mode 
dFomin et al.lll994 ). i.e., by alternately tracking two positions 
at 0.4° offset from the source position. This observing mode 
provides a reliable background estimate for point-like sources 
such as LS 1+61 303. 
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TABLE 1 

Log of the VHE Gamma-ray Observations 



MJD" 


Obs. Time 


Phase Range 


A'(> 300 GeV)" 




(minutes) 






54348. 148 ±0.087 


203.0 


0.470-0.476 


1.1±1.7 


54349. 157 ±0.080 


210.1 


0.508-0.514 


-1.9 ±2.0 


54350. 136 ±0.094 


214.5 


0.544-0.551 


0.2±2.1 


5435 1.1 60 ±0.080 


219.8 


0.584-0.590 


-0.1 ±2.2 


54352.161 ±0.081 


221.1 


0.621-0.627 


15.7 ±2.4 


54353.163 ±0.082 


224.2 


0.659-0.665 


7.7 ±2.1 


54354. 166 ±0.078 


212.5 


0.697-0.703 


2.8 ±2.2 


54355.153 ±0.091 


172.4 


0.734-0.741 


0.3 ±2.3 


54356. 139 ±0.080 


148.7 


0.771-0.777 


0.1 ±2.5 


54357. 153 ±0.092 


178.1 


0.809-0.816 


6.4 ±2.3 


54358. 155±0.091 


179.3 


0.847-0.854 


8.1 ±2.4 


54359. 154 ±0.092 


184.4 


0.885-0.892 


2.9 ±2.4 


54360.154±0.091 


177.3 


0.923-0.929 


5.7 ±2.5 


54361. 154 ±0.094 


183.1 


0.960-0.967 


5.3 ±2.4 


54362.153 ±0.094 


188.7 


0.998-0.005 


-0.6 ±2.4 


54363. 156 ±0.093 


138.6 


0.036-0.043 


7.5 ±2.8 


54364. 149 ±0.099 


195.7 


0.073-0.081 


3.8± 1.9 


Note. — " The central MJD and half-span of each observation is quoted. 



* Uncertainties are given at 68% confidence level. 

The data analysis was performed using t he standard 
MAGIC analysis and reconstruction software (' Albert et alj 
|2008c). The quality of the data was checked and data taken 
with anomalous event rates (very low atmospheric transmis- 
sion, car light flashes, etc.) were rejecte d following the stan - 
dard procedure, as previously done in [Albert etal] (120091) . 
From the remaining events, image parameters were calculated 
jHillas 1985). In addition, the time parameters described in 
lAliu et alj (2009) were calculated as well. For the 7/hadron 
separation, a multidimensional classification procedure based 
on the image and ti ming parameters wit h the Random For- 
est method was used dAlbert et alj|2008dl) . We optimized the 
signal selection cuts with contemporaneous Crab Nebula data 
taken at the same zenith angle range. The energy of the pri- 
mary 7-ray was reconstructed from the image parameters us- 
ing also a Random Forest method. The differential energy 
spectrum was unfolded taking into account the full instrumen- 
tal energy resolution (Albert et al. 2007). We estimate the sys- 
tematic uncertainty to be about 30% for the derived integral 
flux values and ±0.2 for the obtained photon index. For more 
details on the systematic uncertainties present in the MAGIC 
data, see , Albert et al.. (2008c) . 

2.2. X-rays 

We observed LS I h-61 303 with XMM-Newton during 
seven runs from 2007 September 4 to 11, lasting from 12 
to 18 ks (see Table |2]i, amounting to a total observation time 
of 104.3 ks. The EPIC pn detector used the Large Window 
Mode, while the EPIC MOS detectors used the Small Win- 
dow Mode. All detectors used the medium thickness optical 
blocking filter. Data were processed using version 8.0.0 of 
the XMM-Newton Science Analysis Software (SAS). Known 
hot or flickering pixels were removed using the standard SAS 
tasks. Further cleaning was necessary to remove from the data 
set periods of high background corresponding to soft proton 
flares, reducing the net good exposure durations to 67.0 and 
92.6 ks for the pn and MOS detectors, respectively. 

Cleaned pn and MOS event files were extracted for spec- 
tral analysis. Source spectra were extracted from a ^70" ra- 
dius circle centered on the source (point-spread function of 
15") while background spectra were taken from a number of 
source-free circles with ~150" radius (three for the pn de- 



tector, four for MOSl, and five for MOS2). The ex tracted 
spectra were analyzed with XSpec v 12. 3.1 dArnaudl 11996). 
The spectra were binned so that a minimum of 20 counts per 
bin were present and energy resolution was not oversampled 
by more than a factor 3. An absorbed power-law function 
(wabs*powerlaw) yielded satisfactory fits for all observa- 
tions. De-absorbed fluxes in the 0.3-10 keV range were com- 
puted from the spectral fits. 

Additional observations of 2-5 ks each, consisting of sev- 
eral pointings of 0.2-1 .0 ks, were obtained with the Swift/XRT 
from 2007 September 1 1 to 22. The total observation time 
was 28.5 ks. The Swift data were processed using the 
FTOOLS task xrtpipeline (ver 0.12.1 under HEASoft 
6.6). The spectral analysis procedures were the same as those 
used for the XMM-Newton data. Since the hydrogen col- 
umn density was poorly constrained, with typical values of 
(0.5 ±0.2) X 10^2 cm-2, we fixed it to 0.5 x lO^^ cm'^, a typi- 
cal value for LS I H-61 303 also found in the XMM-Newton fits 
and close to the average value for the Swift fits. 

To search for short-term X-ray variability, we also extracted 
0.3-10 keV background-subtracted light curves for each ob- 
servation, binned to 100 s for XMM-Newton and at half-time 
of the sparse 0.2-1.0 ks Swift pointings within each observa- 
tion. For XMM-Newton, we considered the sum of the count 
rate in the three detectors. From these light curves, we com- 
puted the degree of variability as the standard deviation with 
respect to the mean of the count rate divided by this mean, 
and the significance of this variability computed from the 
probability of the count rate being constant. We also com- 
puted hardness ratios as the fraction between the count rates 
above and below 2 keV. 

3. RESULTS 
3.1. VHE Gamma Rays 

The measured fluxes of LS I H-61 303 at energies above 
300 GeV are listed in Table [T] and the corresponding light 
curve is shown in Figure[T](top). The periodically modulated 
peak emission ( Alber t et al. 2009) is prominently seen as the 
highest flux value of these observations at phase 0.62 (fol- 
lowed by a smaller flux at phase 0.66). In addition to this 
established feature we measured, as in 2006 December, sig- 
nificant flux on a nightly basis at phase 0.85. We note that 
there is significant flux in the phase range 0.8-1.0 at the level 
of (5.2 ± 1.0) X 10"'^ cm"^ s"', which is compatible with the 
2a up per hmit we obtained for the shorter 2006 observations 
(Albe rt et al.ll2009l) . 

We also determined the differential energy spectrum from 
the ~10 hr of observations conducted in the orbital phase 
range 0.6-0.7. A power-law fit yields: 

^A^ _ (2.0±0.3stat±0.6syst) X 10"'^ 
dE cm^ s TeV 

/ J. N -2.7±().3„„±().2,y,, 



compatible within error s with our previous measurements 
(lAlbertet al.ll2006ll2009h . 

3.2. X-rays 

We summarize in Table |2] the parameters of the spectral 
fits and variability obtained for both the XMM-Newton and 
Swift/XRT data sets. All X-ray fluxes quoted hereafter are 
de-absorbed. We note that there is no significant hardness 
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TABLE 2 

Log of the XMM-Newton (TOP) AND Swift (BOTTOM) X-RAY OBSERVATIONS 



Exposure Time 



MJD" 


Total 


pn 


MOS 


Phase Range 


f(0.3-10keV)* 


r 


A'h( 


wabs) 


X-/d.o.f. 


Variab./Signif. 




(ks) 


(ks) 


(ks) 




(10"'2 ergcm"2 g-l) 




cm--) 




(%)/((t) 


54347.202 ±0.072 


12.4 


10.8 


12.4 


0.434^0.440 


13.8 ±0.2 (±2.4) 


1.87 ±0.03 


0.517 ±0.014 


446.2/382 


17.1/22.5 


54349. 140 ±0.076 


13.1 


5.6 


9.1 


0.507-0.513 


12.4±0.2 (±2.1) 


1.66 ±0.03 


0.504 ±0.019 


379.8/352 


16.8/15.1 


54350. 199 ±0.106 


18.3 


10.5 


17.1 


0.546-0.554 


13.3±0.2 (±1.6) 


1.66 ±0.02 


0.514±0.013 


536.5/422 


11.8/13.8 


5435 1.160 ±0.099 


17.1 


10.9 


12.1 


0.583-0.590 


13.4±0.2 (±1.6) 


1.68 ±0.03 


0.525 ±0.015 


431.6/400 


12.1/12.8 


54352. 146 ±0.085 


14.7 


5.8 


13.2 


0.621-0.627 


22.9±0.2 (±1.6) 


1.54 ±0.02 


0.538±0.012 


496.7/447 


6.8/ 5.9 


54353. 167 ±0.084 


14.5 


11.9 


14.5 


0.659-0.665 


18.6±0.2 (±1.1) 


1.58 ±0.02 


0.529 ±0.013 


467.5/429 


5.9/ 1.9 


54354. 144 ±0.082 


14.2 


11.5 


14.2 


0.696-0.702 


12.6±0.2 (±1.3) 


1.65 ±0.03 


0.520 ±0.015 


439.5/407 


10.0/ 8.0 


54354.663 ±0.073 


3.3 






0.716-0.722 


10.3 ±0.7 (±2.0) 


1.71±0.16 


0.5 


(fixed) 


5.38/14 


18.4/ 1.1 


54355.670 ±0.070 


2.7 






0.754-0.759 


12.4 ±0.9 (±1.6) 


1.38±0.14 


0.5 


(fixed) 


8.83/15 


10.7/ 0.3 


54356.671 ±0.073 


3.3 






0.792-0.797 


17.6±1.0 (±2.0) 


1.36±0.10 


0.5 


(fixed) 


9.83/26 


9.9/ 0.4 


54357.674 ±0.073 


3.3 






0.830-0.835 


13.9 ±0.9 (±2.3) 


1.49 ±0.14 


0.5 


(fixed) 


13.07/20 


15.3/ 1.2 


54358. 178±0.103 


3.8 






0.848-0.855 


17.7±0.9 (±1.5) 


1.47 ±0.09 


0.5 


(fixed) 


20.46/32 


6.3/ 0.1 


54359.951 ±0.205 


5.3 






0.911-0.926 


19.5 ±0.8 (±2.6) 


1.54 ±0.07 


0.5 


(fixed) 


34.52/50 


12.7/ 0.4 


54362.247 ±0.237 


2.1 






0.996-0.014 


13.7±1.1 (±3.6) 


1.78±0.19 


0.5 


(fixed) 


5.45/11 


25.0/ 1.7 


54363.121 ±0.102 


2.5 






0.034-0.042 


12.3±1.0 (±1.5) 


1.53±0.16 


0.5 


(fixed) 


9.71/13 


9.3/ 0.1 


54365.530 ±0.040 


2.2 






0.127-0.130 


13.3±1.1 (±1.7) 


1.46±0.15 


0.5 


(fixed) 


10.14/12 


10.1/ 0.5 



Note. — " The central MJD and half-span of each observation is quoted. * Fluxes are de-absorbed. The uncertainties in parentheses include the uncertainties 
due to intrinsic X-ray variability during the observations. 



0.4 



0.5 



0.6 



Orbital phase 
0.7 0.8 0.9 



1.0 



1.1 



25 

_ 20 
I y 15 

o t in 
CO o 







25 






20 


> 


w 










Ok 


cm 


15 








cr> 





10 


o 




LL 


O 








5 



+ ff+ \.if ^ 



* ^ 



XMM-Newton 
O— O SwiftlXm 





54345 



54350 54355 54360 

MJD (JD-2400000.5) 



54365 



Fig. 1 . — VHE gamma-ray and X-ray light curves of LS I -h61 303 during 
the multiwavelength campaign of 2007 September Top: flux above 300 GeV 
vs. the observation time in MJD and the orbital phase. The horizontal dashed 
line indicates flux. The vertical arrows mark the times of simultaneous 
VHE gamma-ray and X-ray observations. Bottom: de-absorbed flux in the 
0.3-10 keV energy range for the seven XMM-Newton observations (filled cir- 
cles) and the nine Swift ones (open circles). Error bars correspond to a Icr 
confidence level in all cases. Dotted lines join consecutive data points to help 
following the main trends of the light curves. The sizes of the symbols are 
larger than the time span of individual observations. 

ratio change within each of the observations, in contrast to 
w hat was found in th e XMM-Newton observations reported 
bv lSidoh et alj (120061) . In principle this impHes that, for each 
observation, the flux and corresponding uncertainty obtained 
from the spectral fit is a good estimate of the flux during the 
whole observation. However, moderate count-rate variability 
on timescales of ks is present in most observations, ranging 
from 6% to 17% for the XMM-Newton data and from 9% to 
25% for the Swift data (see Table|2|i. The rms of this variabil- 
ity should be considered in addition to the statistical uncer- 



tainty when providing a flux measurement spanning several 
ks. We converted this count-rate variability into flux variabil- 
ity by m ultiplying the degree of variability defined in Sec- 
tion |22] with the fluxes coming from the spectral fits. Since 
no spectral change is detected within each observation, we 
added this flux variability in quadrature to the spectral fits 
flux errors. This procedure provides the more realistic total 
flux uncertainties quoted in parentheses in Table |2l and used 
hereafter. 

We show in Figure [T] (bottom) the 0.3-10 keV light curve 
of LS I H-61 303. The source displays a steady flux during 
the first four observations and shows a steep increase at phase 
0.62, which is followed by a slower decay up to phase 0.69 
{XMM-Newton) and probably up to phase 0.72 {Swift). The 
behavior is very similar to the one seen at VHE gamma rays, 
although at X-ray energies the baseline has a significant flux. 
Later on there is a significant increase of the X-ray flux up to 
phase 0.8. This high flux is detected with a sparse sampling 
up to phase 0.9, and the source goes back to its baseline flux 
at phase LO. This high X-ray flux between phases 0.8 and 1.0 
occurs when the source is also detected at VHE gamma rays. 

3.3. X-Ray/VHE Gamma-Ray Correlation 

A clear correlation between the X-ray and VHE gamma- 
ray emissions is seen during the outburst, with a simultane- 
ous peak at phase 0.62 (see Figure [T]i. To study the signifi- 
cance of this correlation, we selected the X-ray data sets that 
overlap with MAGIC observations. There are six overlap- 
ping MAGICIXMM-Newton data sets, for which strictly si- 
multaneous observations range from 3.3 to 3.9 hr. For the 
four overlapping MAGlCISwift data sets, the strictly simul- 
taneous observations range from 2.2 to 4.1 hr (although the 
Swift runs have gaps). We plot in Figure |2] the X-ray fluxes 
against the VHE fluxes (from Tables [T] and |2]l for all 10 si- 
multaneous pairs, which are marked with arrows in Figure [T] 
The linear correlation coefficient for the six simultaneous 
MAGIC/XMM-Newton pairs that trace the outburst is r = 0.97. 
For the ten simultaneous pairs, we find r = 0.81 (a \r\ larger 
than that has a probability of about 5 x lO"-' to be produced 
from independent X-ray and VHE fluxes). Minimizing 
we obtain = 7.68 for 8 de grees of freedom and the fol- 
lowing relationship; F(0.3-10 keV)/[10~'^ erg cm~^ s~'] = 
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Fig. 2. — De-absorbed X-ray fluxes as a function of VHE gamma-ray 
fluxes. Only the 10 simultaneous fluxes, marked with arrows in Figure [T] 
have been considered. Error bars correspond to a Icr confidence level in all 
cases. The solid line represents a Unear fit to all data points. 

(l2.2-!^%) + {0Jlt^\l) xN(E >300 GeV)/[lO-^^ cnT^ s"'] 
(non-Gaussian uncertainties). This fit is plotted as a solid line 
in Figure |2] 

However, as can be seen in Figure|2l the flux uncertainties 
are relatively large. This calls for a test of the reliability of the 
correlation strength considering the errors of individual data 
points. To this end, we use the z-Transformed Discrete Cor- 
relation Function (ZDCF), which determines 68% confidence 
level intervals for the c orrelation coefficient from unevenly 
samp led data (see, e.g., l Edelson & Kroliklll988l : lAlexanden 
119971) . The Fisher z- transform of the linear correlation coef- 
ficient is used to estimate the 68% confidence level interval. 
Applying the ZDCF to the X-ray and VHE light curves re- 
ported here we obtain the following uncertainties for the lin- 
ear correlation coefficient: r = 0.81;!;q 2p The sensitivity of the 
MAGIC telescope requires observation times of several hours 
to be able to get few-sigma detections for the VHE fluxes from 
LS I H-61 303. During such long periods, the X-ray fluxes have 
quite large intrinsic variability (see Table |2|. This results in 
quite large uncertainties on both the VHE and X-ray fluxes 
that restrict the ability to determine a possible correlation out 
of the outburst. 

4. DISCUSSION 

We have discovered an X-ray/VHE gamma-ray correlation 
in the gamma-ray binary LS I H-61 303 based on simultaneous 
multiwavelength data obtained with MAGIC, XMM-Newton, 
and Swift during a single orbital cycle. The correlation is 
mainly due to the very similar trends of the detected flux dur- 
ing the outburst around orbital phase 0.62, while the uncer- 
tainties prevent to be sure about the existence of the correla- 
tion outside the outburst. Given the variability in the X-ray 
flux (up to 25% on hour scales), it is necessary to program 
simultaneous observations to perform correlation studies in 
LS I H-61 303. This conclusion has also been reached recently 
by the VERITAS Collaboration, when reporting the lack of X- 
ray/VHE correlation based on contemp oraneous data with a 
VHE sampling that is not dense enough dAcciari et al.ll2009h . 
Although a similar X-ray/VHE correlation has been obtained 
for the VHE gamma-ray emitting X-ray binary LS 5039, this 



result was based on non-simultaneous data acquired y ears 
apart dAharonian et al.ll2005L 120061: iTakahashi et alj|2009l) . 

The VHE emission within a binary system can suffer 
photon-photon absorption via pair creation, mainly with the 
stellar optical/ultraviolet photons. In LS I +61 303, this 
absorption is only expected to be significant tow ard the 
observer for E > 300 GeV just before periastron jPubusI 



l2006btlBednarekl2006l:[s"ierpowska-Bartosik & Torresi2009l) . 
a phase range not explored here. In addition, the quoted X- 
ray fluxes are already de-absorbed (and the hydrogen column 
densities and the associated errors are low). Overall, there 
are no absorption effects to be considered. Therefore, the X- 
rayA'HE correlation we have found for LS 1 +61 303 cannot 
be an artifact due to variable absorption toward the source. 
This indicates that the emission processes at both wavelengths 
occur at the same time and are probably the result of a single 
physical mechanism. In this context, it is reasonable to as- 
sume that the X-ray and VHE emissions are produced by a 
single particle population. 

It is interesting to note that the MAGIC spectrum in the 
0.6-0.7 phase range yields ^ 1 1 x 10"'^ erg cm"^ s"' for E > 
300 GeV, while the X-ray flux is - 19 x 10"'^ erg cm'^ s"'. 
Therefore, the total X-ray flux is approximately twice the 
VHE flux. However, if we subtract an apparent baseline X- 
ray flux of 10 x 10"'^ erg cm"^ s"', the resulting X-ray flux is 
similar to the total VHE flux in the phase range 0.6-0.7. If the 
radiation mechanisms are dominated by a single particle pop- 
ulation, the X-ray/VHE correlation and the smaller/similar 
VHE fluxes favor leptonic models. In hadronic models, the 
X-ray emitting and the VHE photons would come from 
the same protons (for reasonable values of the magnetic field), 
and the luminosity of the radiation should be < 1 /2 that 
of VHE gamma-rays (see Figure 5 of Kelner et al. 2006 for 
reasonable proton energy distributions) unlike it is observed. 
In addition, the inverse Compton (IC) cooling channel is 
less efficient than the synchrotron channel to produce the de- 
tected X-ray emission for reasonable values of the magnetic 
field (see Takahashi et al. 2009 for a similar discussion for 
LS 5039). This clearly suggests that the X rays are the result 
of synchrotron radiation of the same electrons that produce 
VHE emission as a result of IC scattering of optical/ultraviolet 
stellar photons. 

The observed photon indices of the simultaneous X-ray and 
VHE spectra are consistent with one population of electrons 
following a power-law energy distribution with index ^2.1. 
These electrons would produce X-rays via synchrotron and 
VHE photons via IC with an interaction angle < 7r/2. We note 
that an electron index of ^^2.1 is too hard if synchrotron cool- 
ing dominates in the X-ray range, since it implies an injected 
electron index of 1.1. On the other hand, dominant adia- 
batic cooling implies an injection index of 2. 1 , a more reason- 
able v alue (see, for example, the discussion in lTakahashi et al.l 
120091) . Although small changes in the VHE spectrum would 
be expected due to variations in the IC interaction angle or 
the electron index (as seen in X-rays), at present they are not 
detectable due to the large uncertainties of the VHE photon 
index. 

Finally, we note that contemporaneous radio light curves 
obtained with RATAN, VLBA images, and Ha spec- 
troscopy are consistent with previous resu lts (lGregorvll2002l : 
iDhawan etall 120061: IZamanov et all 1 19991) . Details on these 
observations will be reported elsewhere. Therefore, the X- 
ray/VHE correlation occurred when the source was showing 
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a standard behavior in both its outflow (radio) and decretion 
disk (Ha Une). 
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